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We demonstrate the existence of stable three dimensional dissipative localized structures in the 
output of a laser coupled to a distant saturable absorber. These phase invariant light bullets are 
individually addressable and can be envisioned for three dimensional optical information storage. 

An effective theory provides for an intuitive picture and allows to relate their formation to the 
morphogenesis of static auto-solitons and cellular patterns. The complexity incurred by the widely 
different time scales present in the problem as well as by non-local couplings that stem from the 
material degrees of freedom is circumvented by the use of a multiple time-scale analysis. This 
provides a powerful model enabling to tackle effectively the three dimensional case. 


The possibility of using light bullets (LBs), i.e. pulses 
simultaneously confined in the transverse and the prop¬ 
agation directions, has attracted a lot of interest in the 
last twenty years, both for fundamental and practical rea¬ 
sons. Optical confinement is ordinarily envisioned trough 
a mechanism in which a Kerr self-focusing nonlinearity 
compensates the spreading effect of diffraction. Yet in 
this scenario, reminiscent of the conservative soliton the¬ 
ory, conservative LB are unstable and lead to a collapse 
[1] in three dimensions. Other confinement mechanisms 
were envisioned in systems far from equilibrium ande LBs 
were predicted in semiconductor cavities [2, 3]. However, 
it was later shown that the proper consideration of the 
dynamics of the active material leads to collapse [4]. 

In this letter, we establish the existence of stable LBs 
in semiconductor laser cavities. Departing from previous 
approaches, the structuration mechanism exploits the ex¬ 
istence of the active material temporal scales making our 
LB essentially multiple timescale objects. We present an 
intuitive theory that allows to understand the LBs as hy¬ 
brids between the temporal localized structures described 
in [5] and the spatial diffractive autosolitons of [6]. 

Localized structures (LS) have been widely observed 
in nature [7-11] and can be interpreted —in the weak 
dissipative limit— as dissipative solitons [12, 13]. They 
may form when two homogeneous solutions coexist for 
the same values of the parameters [6, 14, 15]. These 
LS correspond to a locking between two opposed fronts 
connecting the two stable states. Another case consists 
in the coexistence of an homogeneous and a modulated 
solution [16] leading to the so-called cellular patterns [17]. 

One separates in optics the systems in which the LS are 
locked to an external injection beam from the ones that 
possess a phase invariance, see [18, 19] for a review. The 
former case leads to the so-called cavity solitons [20, 21] 
observed either in the transverse plane of broad area am¬ 
plifiers [22] or in the temporal output of fibers [23, 24]. 
In the latter case, diffractive autosolitons were predicted 
either in situations where the cavity is composed of a 
gain medium coupled to a saturable absorber [6, 15] or 
to an external diffraction grating [25]. Because these las¬ 
ing LS (LLS) appear in a phase invariant system [26], 
their properties are very different from phase-locked LS, 
leading for instance to optical vortices [27]. 


It was recently shown [5] that phase invariant tem¬ 
poral LLS can evolve from passive mode-locking (PML). 
The PML regime leads to the emission of temporal pulses 
much shorter than the cavity round-trip [28]. It is 
achieved by combining a laser amplifier providing gain 
and a nonlinear loss element, usually a saturable absorber 
(SA). The different dynamical properties of the two ele¬ 
ments create a window for regeneration only around the 
pulse. Although many pulses may coexist within a PML 
cavity, they cannot be addressed individually. It is due to 
the fact that the dynamics cannot be reduced to the sole 
evolution of the field and while the pulses are typically 
short and may seem temporally localized, they may still 
strongly coupled trough the evolution of the gain that 
occurs on a much longer time scale. 

For cavities with a large aspect ratio, as defined by the 
ratio of the gain recovery time and of the cavity round- 
trip time r = r/r^, the PML pulses may become under 
certain conditions individually addressable LLS coexist¬ 
ing with the off solution [5]. In this regime of temporal 
localization, the LLS were found to arise from a bifurca¬ 
tion scenario radically different to the one found in PML 
lasers. While PML occurs as a supercritical Andronov- 
Hopf bifurcation over a continuous wave solution, here 
the LLS are nascent from Saddle Node bifurcations of 
Limit cycles (SNL) and occur below the lasing thresh¬ 
old. As such, they coexist between themselves as well as 
with the zero intensity solution. The coexistence with 
the harmonic PML solution of maximal order hinted at 
a similarity with decomposable cellular patterns [17]. 

In this manuscript we discuss in which conditions the 
transverse profile of these temporal LLS can self-organize 
yielding to a new robust LB formation scenario. We de¬ 
scribe the PML laser using the generic delayed differential 
equation model presented in [29] which can describe both 
the pulsating regimes and the steady solutions. We work 
in the limit of low gain (G) and saturable absorption (Q) 
as to justify a first order approximation to the single pass 
evolution of the field profile. In addition, we assume that 
the transverse section is sufficiently broad for the effect 
of diffraction at each round-trip to be small. This allows 
the lumping and the commuting of the various nonlinear 
elements. In this uniform field limit, the equation for the 
field amplitude reads 
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where 7 is the bandwidth of the spectral filter, = 
transverse Laplacian, n is the fraction of 
the power remaining in the cavity after each round-trip 
and a and (3 are the linewidth enhancement factors of the 
gain and absorber sections, respectively. In Eq. (1), the 
transverse space variables r±_ = (x, y) have been normal¬ 
ized to the diffraction length. As such, the domain size 
representing the dimension of the broad area laser 
becomes a bifurcation parameter: we foresee that with 
Lx ^ 1 , one may only find a uniform spatial state while 
localized patterns may occur when L_\_ ^ 1. 

In addition to the two reversible transverse spatial di¬ 
mensions r±_ (found setting r± —^_l), the delayed 

values of the variables (E^G^Q) render Eq. (1) infinite¬ 
dimensional along the propagation direction. The time 
delay (r) describes the spatial boundary conditions for a 
loop cavity and governs the fundamental repetition rate 
of the PML laser. The carrier equations read 

dtG = rGo-G(T + \Ef^+VgA^G, ( 2 ) 

dtQ = Qo — Q (l + s + 'DqAj_Q, (3) 

with Go fhe pumping rate, T = the gain recovery 
rate, Qo is the value of the unsaturated losses which de¬ 
termines the modulation depth of the SA and s the ratio 
of the saturation energy of the gain and of the S A sections 
while the scaled diffusion coefficients are Vg and Vg. In 
Eqs. (1-3) time has been normalized to the SA recovery 
time. The slow temporal scales associated with G and 
Q prevent their adiabatic elimination which renders the 
dynamics along the propagation (time) axis irreversible, 
notwithstanding the presence of dissipative terms. 

The lasing threshold above which the off solution, de¬ 
fined as (L,G,Q) = (0,Go,Qo), becomes unstable is 
Gth = — 2 + Qo. We study the dynamical sys¬ 

tem defined in Eqs. (1-3) in the newly found regimes of 
LLS [5], where the cavity round-trip time is much longer 
than the gain recovery r = r/r^ ^ 1 and for values of 
the gain below threshold, i.e. Go < Gth- In this regime, 
LLS occur essentially as multiple timescales objects due 
to the widely different dynamics of the light and of the 
matter components. 

We show in Eig. 1 that Eqs. (1-3) can indeed support 
the emission of LBs as temporal LLS that are also con¬ 
fined in the transverse plane. Here the extent of the 
optical component of the LB normalized to the round- 
trip time is typically r^/r ^ 10“^. Due to the very large 
aspect ratio of our system, such LBs would be invisible 
in a spatio-temporal representation for the optical field 
intensity as it would corresponds to an extremely thin 
horizontal segment. It is why we represented in Eig. 1 
the slower material variables G and Q. Starting from the 
off solution, i.e. the solution with N = OLB, each ad- 
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Figure 1. (Color online) Spatio-temporal trace for the gain 
G (x, t) (a,c) and for the absorption Q (x, t) (b,d) in a regime 
where N — ?> LBs are present. The time windows is such 
that each LB appears twice. The inset (b,d) highlight the 
difference between the relaxation timescales of the material 
variables. The slowest recovery variables G defines the tem¬ 
poral extent of the LB. The parameters are 7 = 40, n — 0.8, 
a = /3 = 0 , T = 200 , Go = 0.6678Gt/^, F = 0.04, Qo = 0.3, 
s — 30, Lx = 90 and Vg^q — 0. 


ditional LB was triggered individually by sending a per¬ 
turbation in the form of an optical pulse. Since LBs are 
attractor, the details of the writing pulse are irrelevant 
as long as its energy remains above a certain threshold. 
After a transient of several tens of round-trips, the sit¬ 
uation reaches an equilibrium and an additional LB is 
written. With respect to the existing LBs, a new LB can 
be written either at the same time but at a different lo¬ 
cation, like the two LBs represented in the white inset in 
Fig. 1, or a two different times. In this case, then can 
even partially overlap spatially like the two leftmost LBs 
in Fig. 1 at coordinates (x,t/r) = (-40,-0.75). With¬ 
out this staggering, these two LBs could not coexist at 
the same instant. 

In the case of a single LB, we performed a bifurcation 
diagram as a function of the gain parameter. Our re¬ 
sults are depicted in Fig. 2 where we represent the pulse 
energy. Here one observes that the LBs occurs below 
the lasing threshold. The square root behavior around 
the minimal gain value suggests a scenario based upon a 
SNL bifurcation. For high current, the LB spatial profile 
develops wings that oscillate at a low frequency signaling 
the onset of an Andronov-Hopf bifurcation. We depict in 
Fig. Ib-e) the details of the LB at the lower and highest 
stable currents. 

We stress that the self-organization mechanism at work 
is not based upon a self-focusing Kerr nonlinearity. For 
instance, we obtained the results of Fig. 1 and Fig. 2 
setting the values a = /3 = 0. Exploiting the seminal 
work of New [30] and the fact that the LLS are com- 
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before. In our case where r ^ Tg^ the material variables 
perform a full recovery in-between the emission of a LLS. 

We wrote Eq. (4) to make apparent the link between 
our approach and the work of [31, 32] for the case of static 
auto-solitons in bistable interferometers. In [31, 32], one 
assumes a monomode continuous wave emission along the 
longitudinal propagation direction, which allows, via the 
adiabatic elimination of the material variables, to find an 
effective equation for the transverse profile, yet with a 
different expression for the function /. This would corre¬ 
sponds in Eqs. (1-3) to taking the limit r ^ 0 and setting 
dtG = dtQ = It is striking that in the frame of a tem¬ 
poral train of short pulsating LLS, one finds exactly the 
same effective equation describing the transverse struc¬ 
turation dynamics. Eollowing the method detailed in [32] 
in the case of a single transverse spatial dimension, we in¬ 
troduce in Eq. (4) a phase-amplitude decomposition and 
a spectral parameter uj as 

^4 (x, cr) = p (x) exp {i [0 (x) — cjcr]} , (6) 

and are left searching for static spatial LLS as hetero¬ 
clinic and homoclinic orbits of 

dxP = kp , dxq =-2qk + lR[f (7) 

dxk = -w + - 9 [/ (p)] , (8) 


Figure 2. (Color online). Bifurcation diagram for the spatially 
integrated peak power of the LB intensity (P) as a function 
of the gain Gq. The single LB regime occurs as a SNL bifur¬ 
cation around Gsn — 0.668Gth and looses its stability via an 
Andronov-Hopf bifurcation at Gh ~ 0.694Gt/j. In both cases, 
the spatio-temporal profile of the field is represented in pan¬ 
els b,d) and c,e) with Gq — Gsn and Gq — Gh, respectively. 
The others parameters are like in Fig 1. 


posed of different variables evolving over widely different 
timescales, we find an effective equation for the trans¬ 
verse morphogenesis problem in presence of diffraction. 
We assume that the field reads E (r_L, t) = A (r_L, t)p (t) 
with p{t t) = p{t) a short normalized temporal pulse 
train and A (r_L, t) a slowly evolving amplitude. Inserting 
this Ansatz in Eqs. (1-3), defining a = yt, we find that 
the equation governing the dynamics of the transverse 
profile reads 


d,A = iA^A + Af(\Afy (4) 
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while the expression of the nonlinear function / reads 
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To obtain Eqs. (4,5), we used the fact that, during the 
emission of a LLS, the stimulated terms are dominant 
in Eqs. (2-3) which allows connecting the values of the 
material variables after the pulse impinging to the values 


Figure 3. Bifurcation diagram of the reduced model Eq. (4) 
as a function of the gain parameter Go. The panels a,b) rep¬ 
resent the evolution of the spectral parameter uj while the 
spatially integrated intensity S is depicted in panels c,d). The 
strong multistability of solutions is depicted in panel e). Usu¬ 
ally, only the transverse LS of higher intensity is stable. The 
main solution is stable in the vicinity of the Saddle-Node bi¬ 
furcation that occurs at Gsn = O.GOGGth and losses its sta¬ 
bility via and Andronov-Hopf bifurcation at Gh = 0.63Gt/i. 
We represent the oscillations of the transverse profile slightly 
above Go — Gh in panel f). Parameters as in Fig. 1. 
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Figure 4. Bifurcation diagram as a function of Gq. We repre¬ 
sent the drift velocity of the LB (a), their energy (b) as well 
as the field (c), gain (d) and absorption (e) spatio-temporal 
profiles at Go = 0.674Gth, indicated in panels a,b) by a circle. 


where we defined k = p~^dxp and q = Upon con¬ 

tinuation as a function of a parameter one is able to 
reconstruct a full bifurcation diagram and our results 
are shown in Fig. 3. Here, the spiral shape of the sin¬ 
gle LS branch is reminiscent of the results of [32] in the 
case of static auto-solitons. The stability of the solution 
branch was found via the reconstruction of the linear evo¬ 
lution operator along the lines of [33, 34]. As detailed in 
[32], Eq. (4) exhibits three zero eigenvalues related to 
the phase, translational and Galilean invariance. This 
reduced model bifurcation diagram is in good qualitative 
agreement with one of the full model. The extent of the 
stable LB region is well reproduced as AG ^ 0.03Gt/i. 
The rigid shift to the left in the reduced model can be 
ascribed to the neglect of the spectral filtering loss term 
inherent to New’s approach. Finally, the case of radially 
symmetric 2D LBs can be studied similarly. 

However, the aforementioned reduced model does not 
allow access to the temporal pulsewidth. One can trans¬ 
form Eqs. (1-3) into a more tractable PDE using a multi¬ 
ple time scale analysis. In the limit of a broadband filter 
7 ^ oc, we reduce Eqs. (1-3) to the Hans equation 


d,G = VGo-G{T + \E\^'^+Vg/\^G, (10) 

dzQ = Qo — <3 ^1 + s |-B|^) + 'DqA^G, (11) 

with T a slow time scale. The Eqs. (9-11) clarify the com¬ 
posite nature of the LBs as one notices the presence of 
diffusion along the longitudinal dimension and diffraction 
in the transverse dimensions. Although Eqs. (9-11) may 
seems a complication with respect to Eqs. (1-3) since the 
system is now 4D, it allows to cut the domain along the 
2 )-axis to a box of the size of the optical pulse and neglect 
the long tail of the gain recovery. The reason for doing so 
is that the gain material losses entirely its memory at the 
next round-trip, so that one can always set the bound¬ 
ary condition (G, Q) {z = 0,r_L,T) = (Go,Qo)- One can 
appreciate in Eig. 4 the results of this approach and no¬ 
tice the excellent agreement with the results of the full 
model. Einally, this approach can readily be extended to 
the study of 3D LB as depicted in Eig. 5. 

In conclusion, we have shown that broad area PML 
lasers with a large temporal aspect-ratio can display ad¬ 
dressable Lasing Light bullets. We presented an intu¬ 
itive theory that allows to understand these multiple time 
scales object as hybrids between the temporal localized 
structures [5] and spatial diffractive auto-solitons [6]. A 
multiple time scale analysis allows reducing the size of 
the longitudinal domain by several orders of magnitude 
rendering the three dimensional analysis feasible. 
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Figure 5. Isosurface representation of a three dimensional LB 
at which the Intensity (a), the real (b) and the imaginary 
parts (c) are 1 % of the maximal value. The multiple curves 
in (c) indicate that the LB are spatially oscillating. 
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